Thermalization in open classical systems with finite heat baths 
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We discuss thermalization of a test particle schematized as a harmonic oscillator and coupled to 
a Boltzmann heat bath of finite size and with a finite bandwidth for the frequencies of its particles. 
We find that complete thermalization only occurs when the test particle frequency is within a 
certain range of the bath particle frequencies, and for a certain range of mass ratios between the 
test particle and the bath particles. These results have implications for the study of classical and 
quantum behaviour of high-frequency nanomechanical resonators. 
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Since the pioneering studies on Brownian motion and 
its interpretation at the molecular level [H, the study 
of open systems has been a crucial area in both classi- 
cal and quantum statistical mechanics. The theory of 
open systems allows one to bridge the closed, pure-state 
dynamics of a particle with the open, mixed-state dy- 
namics in the presence of coupling to an environment. 
The deterministic dynamics of the particle considered 
as a closed system is replaced in the quantum regime 
by a stochastic Schrodinger equation, corresponding to 
a stochastic Newton equation in the classical limit. In 
the classical case, the environment can be represented 
as a heat reservoir consisting of a large set of harmonic 
oscillators with Boltzmann energy distribution, all sym- 
metrically coupled to a specific particle, which hereafter 
we will refer to as the test particle. In the limit of an 
infinite number of oscillators, and by properly choosing 
the density of states and the initial conditions, one can 
derive a stochastic equation which leads to solutions with 
Boltzmann energy distribution for the test particle, i.e. 
thermalization at the reservoir temperature. 

To the best of our knowledge, however, there have been 
no investigations of the conditions under which thermal- 
ization of a test particle occurs in the presence of a re- 
alistic heat bath with limited resources; namely a finite 
number of degrees of freedom and a finite bandwidth. 
The closest discussion in the literature regards thermal- 
ization of a test particle with an independent oscillator 
heat bath, but with oscillators distributed in a generic 
and unbound frequency range Q , the study of interact- 
ing Fermi systems with finite numbers of particles [1] , and 
the celebrated Fermi-Pasta- Ulam problem, where nonlin- 
ear interactions are crucial [J . To consider a finite band- 
width for a finite number of oscillators making a heat 
bath is significant because in many contexts, in partic- 
ular for systems of interest in mesoscopic physics and 
nanotechnology, the size of the environment is small and 
may no longer justify the thermodynamic limit. This 
leads naturally to an infrared cut-off for the density of 
states of the bath, as the latter cannot support wave- 
lengths much larger than its spatial extent. Furthermore, 



the finite amount of energy in any realistic environment 
demands an ultraviolet cut-off in its density of states. In 
this paper, we describe results on the thermalization of 
a test particle in the presence of a heat bath comprising 
a finite number of harmonic oscillators with frequencies 
distributed within a finite bandwidth. The main result 
of our analysis is the emergence of diverse scenarios for 
the statistical distribution of the energy of the test par- 
ticle, such as thermalization at the heat bath tempera- 
ture, Boltzmann distribution at a smaller temperature 
than the heat bath, or no thermalization, depending on 
the interplay between the intrinsic frequency of the test 
particle and those of the heat bath. While our study fo- 
cuses on the classical case, it points to a broad range of 
implications in the physics of nanomechanical structures, 
whose displacements may be monitored at or beyond the 
standard quantum limit 0,0]. Furthermore, some of the 
results reported here may find application wherever an 
open system approach is necessary, encompassing sys- 
tems as diverse as classical plasmas and self-gravitating 
objects, ultracold trapped atoms and quasiparticles in 
nanostructures. 

The starting point for our analysis is the Hamiltonian 
of the test particle plus environment 0, S, 0] : 
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Here the test particle, with generalized coordinates 
{Q,P), is modeled as a harmonic oscillator of mass M 
and angular frequency il. The n-th particle in the reser- 
voir, a harmonic oscillator of mass m and angular fre- 
quency LUn, is characterized by generalized coordinates 
ilniPn)- To avoid unnecessary complications at this 
stage, the systems are assumed to be one-dimensional. 
The coupling between the test particle and the parti- 
cles in the reservoir is chosen to be translationally in- 
variant, which eliminates the appearance of renormaliza- 
tion terms l3, U]- By writing the Hamilton equations 
and solving for the particles of the heat bath, a general- 
ized Langevin equation for the test particle is obtained 
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dsT{t - s)Q(s) + n^Q(t) = n(t), (2) 



where the kernel of the dissipative term (non local in 
time) and the fluctuation force term are respectively: 
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In the limit of inflnite N extending over the contin- 
uum frequency range ujn S [0, +oo). with density of states 
dN/dco oc and if the initial conditions for the bath 

particles {qn{to),Pn{to)) are chosen to realize a Boltz- 
mann energy distribution, the test particle is described 
by a Langevin equation (for an extensive discussion see 
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dQ{t) = P{t)dt/M, 

dP{t) = -[7P(t) + iim2Q(t)]dt, 
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where 7 = {'nmuP' l2M){dN / du) is the dissipation factor 
(which becomes independent of u) for a density of states 
(X 1/w^) and 5 = y/2m^kBT. The stationary solution of 
Eqs. (5)-(6) corresponds to the thermalization of the test 
particle energy with the same Boltzmann distribution as 
the bath. Such a situation is highly idealized, since it 
corresponds to an infinite number of oscillators with all 
possible frequencies, and these conditions are not always 
met in real physical systems. 

To study the conditions under which thermalization 
occurs in a more realistic setting, we have investigated 
the dynamics of a reservoir characterized by a finite 
frequency spectrum a;„ € [wirjWuv], a finite number 
of oscillators N, a finite mass ratio m/M between the 
masses of the bath particles and the test particle, for a 
uniform density of states, obtained by a random choice 
of the bath particle's frequencies, such that dN/dtu cx 
9{lo - Wir) + 9{u;uv - lo) - 1. 

We do this by generating explicit numerical values for 
the parameters which describe each bath oscillator, solv- 
ing the equations of motion from Eq. (1) for the motion 
of the test particle, and sampling the test particle energy 
over time. The initial conditions are chosen, for the bath 
oscillators, to be randomly distributed in phase space in 
such a way as to generate a Boltzmann distribution for 
their energy, and consequently a well defined tempera- 
ture, while the test particle starts with zero total energy. 
The system of iV-f 1 oscillators can then be solved exactly 



after numerical diagonalization of the matrix describing 
the Hamiltonian evolution. Since the test particle is a 
single system, no meaningful definition of its tempera- 
ture may be given in terms of the energy distribution 
of a statistical ensemble at a given time. An effective 
temperature can still be defined by sampling the energy 
of the test particle during its time evolution, collecting 
an energy distribution over time, and fitting the energy 
distribution with a Boltzmann distribution law. This 
operative definition does, however, preclude observation 
of transient behaviour, at least on timescales smaller or 
comparable to the time interval over which the energy 
distribution is collected, in the test particle temperature. 
Indeed, we find that the first effective temperature we are 
able to define, when enough sample points have been col- 
lected to limit the statistical error, is already the equilib- 
rium temperature. This precludes in principle the study, 
for instance, of the dynamics of nonequilibrium station- 
ary states, nonequilibrium phenomena such as aging in 
glassy systems, or the approach to thermal equilibrium. 
Also, to avoid the effect of possible recurrences and spu- 
rious periodicities we have randomly chosen the sampling 
times for the measurement of the test particle energy. 

The energy distribution of the test particle is shown 
in Fig. 1 for various values of the bandwidth Aw = 
ujuy — oJiR of the frequency spectrum of the reservoir. 
For a test particle interacting with a finite number of 
oscillators all degenerate in frequency (a;„ = ujr), we 
do not observe thermalization, as the resulting energy 
distribution is quite far from being Boltzmann in nature. 
Indeed, if the initial conditions for the bath oscillators 
are distributed symmetrically in phase space (such that 
Enlni-to) - 0, EnPni-to) - 0), Eqs. (3)-(4), and (2) 
may be approximated respectively as: 
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Q{t) + ^uj% Jl Q{s) sm[LUB,{t - s)]ds + tt^Q{t) = 0,(9) 

where Cl = {1 + ^y/^fl and ( ^ Nm/M. 

The dynamics corresponding to Eq. (9) resembles that 
of multiple oscillators experiencing beating phenomena, 
with energy periodically transferred between the test par- 
ticle and the heat bath with a harmonic law, as iStp (t) = 
i;o{l+C/2[l-cos(m/2)]}, where is the bath's initial 
energy. The corresponding probability density P{E) for 
the test particle to have energy between E and E + dE is 
P{E) cx (dE/dt)-^ cx {{E-EQ)[EoilH)-E]}~^/^. This 
behaviour is shown in Fig. la; evidently a monochro- 
matic reservoir is too deterministic to allow scrambling 
of the exchanged energy, and we see here that the test 
particle does not thermalize. With even a modest band- 
width for the reservoir the coherent driving of the test 
particle does not take place at the same level, as shown 
in the top-right and bottom-left plots for increasing val- 
ues of the bandwidth, and above a given value the energy 



3 




Energy 



Energy 



FIG. 1: Dependence of the energy distribution of the test 
particle on the frequency bandwidth of the bath osciUators. 
The number of times A^tp the test particle is found to be with 
energy between E and E + AE (with AE = 0.5 the energy 
binning interval) is plot versus the energy of the test parti- 
cle (expressed in arbitrary units). The heat bath consists of 
400 harmonic oscillators with frequencies centered around the 
test particle frequency, f2 = 60 in arbitrary units, and with 
bandwidth equal to Aco =0 (resonant bath, (a)), Aoj =10 (b), 
Alj =40 (c), and Au; =80 (d). The test particle energy, ini- 
tially assigned to be zero, is sampled at random times, on av- 
erage every 0.1 times its intrinsic period of oscillation. Apart 
from the degenerate bath case [l^ , we have observed that the 
test particle approaches its stationary energy state regardless 
of its initial conditions provided that these correspond to an 
initial energy small enough. Measurements on the test par- 
ticle are taken after waiting a time long enough to yield sta- 
tionary distributions, with a total of 10* measurements. The 
mass ratio between the bath particles and the test particle is 
m/M = 10"^ 



distribution eventually approaches the Boltzmann distri- 
bution, as shown in Fig. Id. 

These results are obtained by keeping the spectrum of 
the reservoir centered around the frequency of the test 
particle. An asymmetric situation, in which the oscilla- 
tion frequency of the test particle is not centered within 
the heat bath bandwidth, is depicted in Fig. 2. 

One can obtain situations with no thermalization, ther- 
malization at a temperature different from the one of the 
reservoir, or complete thermalization, depending upon 
the relative location of the test particle frequency with 
respect to the frequency spectrum of the reservoir. The 
reservoir, whose energy distribution is depicted in the 
upper-left panel, is chosen with a bandwidth equal to 
the one in Fig. Id. In the other five plots in Fig. 2, 
the energy distribution of the test particle is shown for 
increasing values of its oscillation frequency. In partic- 
ular, at < wiR (panels (b) and (c)), the test particle 
does not thermalize, and its energy distribution shows 
a peak at a finite energy. In fact, in the limit f2 — > 
(free test particle) and for Auj (resonant bath) the 
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FIG. 2: Dynamics of thermalization. From top-left to bottom 
right: energy distribution for the heat bath composed of 400 
oscillators (a), and energy distribution of the test particle, 
after 10** random samplings, with a frequency Q = (iJiR/20 
(b), Q = oJiR (c), n = (oiiR + ujuv)/2 (d), = cjuv (e), and 
Q = 20aJuv (f). The straight lines are best fits with Boltz- 
mann distributions. Since, as visible in (b) and (c), the energy 
distribution differs significantly from the Boltzmann distribu- 
tion at small n, the corresponding quoted temperatures from 
the Boltzmann fitting have to be interpreted as effective tem- 
peratures, measuring the average energy of the test particle 
rather than the slope of a Boltzmann distribution. 



equations of motion correspond to a driving force and an 
unlimited increase of energy in time. The presence of a 
finite Auj introduces dephasing among the oscillators and 
consequently a finite and stationary energy peak propor- 
tional to ^. In panel (d), the particle reaches a Boltzmann 
distribution with a temperature close to the one of the 
reservoir Ttp ~ Tics, whereas in panels (e) and (f) the 
Boltzmann distributions of the test particle correspond 
to significantly lower temperatures. 

By repeating the procedure for different values of the 
angular frequency of the test particle, we find three dif- 
ferent regimes, as shown in Fig. [3l At low angular 
frequencies << ljjr, the test particle approaches an 
equilibrium state of non-Boltzmann nature correspond- 
ing to an effective temperature significantly lower than 
the reservoir temperature, with a flat dependence upon 
the frequency. The lack of thermalization in this situa- 
tion resembles the situation of the Brownian motion of a 
free particle (obtainable in our case as the limit ^ 0) 
coupled to a bath with super-Ohmic density of states 
In the intermediate regime wir < fl < wuVj the test 
particle generally approaches complete thermalization at 
Ttp ~ Tros. In the region where CI » touy, the test parti- 
cle has a Boltzmann energy distribution at a temperature 
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FIG. 3: EfRciency of thermalization versus test particle fre- 
quency. Shown is the ratio between the temperature of the 
test particle and the temperature of the heat bath versus the 
angular frequency of the test particle (in arbitrary units), for 
different mass ratios m/M. The heat bath is composed of 400 
oscillators with uniform density of states between tJiR=20 and 
ci^uv^lOO, as indicated by the brown rectangle in the lower 
part of the plot. 
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FIG. 4: Dependence of the efficiency of thermalization upon 
the particle masses. Shown is the ratio between the tempera- 
ture of the test particle and the temperature of the heat bath 
versus the ratios between the mass of the test particle and 
the mass of the oscillators in the heat bath. Different plots 
correspond to different frequencies of the test particle, from 
Q, = O.lwiR to SI = lOwuv, as described in the legend. 



Ttp << T^es "which is frequency-dependent. This may be 
understood by considering the fact that a finite ultravi- 
olet cut-off for the heat bath imphes non-Markovian be- 
haviour. Then a high-frequency test particle should be 
dynamically decoupled from the heat bath in analogy to 
the mechanism introduced in 17, 1^ for open quantum 



systems. Notice that, in the case of m/M = 0.1 in Fig. 
3 (and, less manifestly, for the curve at m/M = 10~^), 
the region of nearly complete thermalization is shifted 
to"wards significantly larger frequencies "with respect to 
the interval [wir, a;u"v]- The obtained data are stable 
"with respect to the choice of the initial time at "which 
the sampling of the test particle energy is performed; 
even "waiting more than four orders of magnitude does 
not change the thermalization pattern. 

This behaviour may be qualitatively understood by 
thinking of the case of a degenerate heat bath. In this 
extreme situation, the particle is either in resonance "with 
the bath oscillators, trading energy in an efficient "way, or 
other"wise al"ways detuned, limiting the energy exchange 
both in amplitude and speed. Since the effective fre- 
quency of oscillation of the test particle in Eq. (9) is 
renormalized as 17 = (1-1- ^)^/^r2, if ^ is significantly 
larger than unity "we expect the resonance to occur at 
angular frequencies higher than the intrinsic angular fre- 
quency of the test particle. In the cases of m/M = 10~^ 
and m/M ~ 10^^ we expect, respectively, Vl = \/5fi 
and f) = -\/4T0, and the shifts of the thermalization re- 
gions are in agreement "with this scaling. This has as a 
further implication a non monotonic dependence of the 
thermalization upon the mass ratio m/M, "which cannot 
be trivially addressed in a perturbative approach [l9| . 



At constant f2, we expect that by increasing the mass 
ratio, the parameter ^ becomes so large as to detune the 
effective response of the test particle, making its ther- 
malization to the heat bath less efficient. 

This feature is clearly evidenced in Fig. [4l where the 
mass ratio dependence is shown for four frequencies in 
diiferent spectral regions. For small values of fi, ther- 
malization is initially proportional to m/M, then de- 
creases when this parameter is in the 10^^ range. For 
Q. » wuv, the behaviour of Ttp/Tj-os is well described 
by a direct proportionality law, i.e. Tip/Treg — m/M. 
We have also repeated the simulations for various exam- 
ples of the density of states of the particles in the heat 
bath, with very little effect on the thermalization features 
apart from the peak region, as will be reported elsewhere 



In conclusion, we have studied thermalization of a test 
particle in contact with a finite reservoir, showing that 
this occurs only for specific relationships between the fre- 
quency and the mass of the test particle and the cor- 
responding quantities of the particles constituting the 
reservoir. These results may be understood as the in- 
termediate case between the two extreme situations of a 
reservoir with an infinite bandwidth and a reservoir with 
a degenerate spectrum of frequencies. The analysis can 
be extended to generic, nonlinear couplings between the 
test particle and the reservoir with numerical techniques. 

The interpretation of two recent experiments in 
nanomechanics may be reconsidered in light of the results 
obtained here. In particular, the observation of discontin- 
uous features in the mechanical transfer function of high- 
frequency nanoresonators operated at low temperatures 
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reported in 2l|, could be attributed to non-equilibrium 
features of the test particle, in this case a high-frequency 
mode of a nanomechanical structure. Our conjecture is 
that for nanomechanical structures the infrared cut-off 
frequency wir (ultraviolet cut-off frequency) is signifi- 
cantly larger (smaller) than in the bulk. The infrared 
cut-off is determined by the lowest frequency which sus- 
tains stationary vibrations in the structure, and for typ- 
ical sound speeds {vs — 10^ — 10"* m/s) and typical sizes 
involved in high frequency nanoresonators (100 nm-l/im) 
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2l| . we can easily have wir, ~ 1-10 GHz, realiz- 
ing a situation where < wir. Likewise, the Debye 
frequency could differ significantly with respect to the 
bulk situation due to the effective lower dimensionality 
of the nanostructure. The effective temperature felt by 
the nanoresonator could then be quite different from the 
one expected by measuring the external bath to which it 
is coupled. This should be confirmed by both dedicated 
experiments on thermal properties, such as specific heat 
and heat conductivity, ab initio calculations of the den- 
sity of states, as well as repetitions of the experiments 
with higher frequency resonators and modified size of the 
structures, possibly leading to a new route to investigate 
non-equilibrium quantum statistical mechanics 22 1. 

Also, in [i^, a phenomenon of cooling has been in- 
terpreted as due to quantum back-action of the read-out 
system, even though the system was well above the stan- 
dard quantum limit. A more economical interpretation 
is instead available by imagining the effective tempera- 
ture of the nanoresonator as resulting from the competi- 
tion between two effective classical heat baths at differ- 
ent temperatures. This is in line with the phenomenon 
of cold damping introduced eight decades a go [24l| . and 
demonstrated for macroscopic resonators in [25j . More- 
over, the data show intriguing anomalies, since the effec- 
tive noise temperature and the noise spectral densities 
reported give values 10^ and 15 times above the stan- 
dard quantum limit, respectively [2^. This large dis- 
crepancy could be due to a partial thermalization of the 
nanoresonator in the presence of two competing baths. 
From our perspective, the study of the thermalization 
of a particle in simultaneous interaction with two heat 
baths should shed light on this phenomenon already at 
the classical level, especially focusing on the energy dis- 
tribution which, in the presence of finite resources, is not 
necessarily of Boltzmann nature. 
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